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CLOSED FORM SATELLITE TRACKING DATA CORRECTIONS 
FOR AN ARBITRARY TROPOSPHERIC PROFILE 

John W. Marini 
ABSTRACT 

The formulas commonly used to correct satellite tracking data for the effects of 
atmospheric refraction are based, for the sake of mathematical convenienc , on 
simple model atmospheric profiles that are not physically realistic. A method 
is given here for deriving correction formulas that do not suffer from this limi- 
tation, and can be used in precision orbit calculations. The refractivity of the 
atmosphere is assumed to liave spherical symmetry, but may have any given 
vertical profile. The method was tested for numerical accuracy by application 
to the simple exponential profile, and the corrections calculated agreed closely 
at all elevation angles from 0 to 90 degrees with those obtained by double- 
precision ray-tracing. The error was in all cases less than 1%, and less than 
1/3% above 1 degree elevation. 

The method has been used to obtain improved correction formulas, and these 
will be published in a separate report. 
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CORliEC'l'ING SA'I’ELLI I E 'rUACKlNG DATA FOR A 


SrniERICALLY-SYAIME'rRlC ATMOSPHERE 


1. INTRODLTCTION 


1.1 Background 


In the radio tracking of a satellite by a ground station, measurements are made 
of the satellite elevation, range, and range- rate. The elevation can be determined 
by a measurement of the angle- of- arrival of a radio wave from the satellite. 

The range is obtained from a measurement of the delay of signals propagated 
between the ground station and the satellite, while the range- rate is obtained by 
counting cycles of the Doppler-sliifted received signals over a prescribed s o riod 
of time, or, alternatively, by noting the time interval required to count a pre- 
scribed number of cycles. 


The passage of the radio waves through the atmosphere of the earth introduces 
into these measurements errors that may require correction. Both tropospheric 
and ionospheric effects are present. However, if the radio frequency used is 
sufficiently high, the ionospheric correction may either be neglected or super- 
posed on the tropospheric correction 1 


1.2 Related Work 

A formula for the bending of a radio ray in an exponential atmosphere as a func- 
tion of the angle-of-arrival has been given by Thayer ' 2 . Freeman 3 has 
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obtained, for an exponential atmosphere, a range~error correction as a function 
of the elevation angle of the satellite. The correction is of the first order in the 
surface refractivity. Reichley I 4 J has obtained both second order elevaticn- 
error and range-ei ror correction formulas as functions of tho elevatio.i angle. 

His formulas can be applied to profiles other than exponential. HopReld I 5 ] , by 
neglecting (Xith curvature, has obtained range and range- rate corrections for 
the two-quartic tropospheric refractivity profile. Rowlandson and Moldt 1 6 I 
have derived closed-form corrections for an exponential atmosphere. 

1.3 Satellite Tracking Data and Correction Formulas 

The satellites tracked are ordinarly above the region, which extends from the 
ground to about 70 km in altitude, where most of the radio- ray bending caused 
by the nonionized atmosphere of the earth takes place. 

In a typical pass of satellite over a ground station, the satellite might be under 
observation for a number of minutes. During this period of time, angle-of- arrival 
range and range- rate measurements are taken periodically at a rate, perhaps, 
of one set of data per second [7 1. Consequently, a large volume of data is gener- 
ated, which must be processed automatically by electronic computer. 

To be usablo in a practical sense, the equations employed to correct this data 
for the effects of the atmosphere should require a minimum of computer time. 








It is not easy, h«nvever, ti» dcjive < on’tvtion cqiiationH tluit, on the one liand, 


accurately represent tlic atmospheric ijintjc ! uSHiiined uiid, Oii the other hand, 
are suitable lor use in tlie i>rocessine el large cjuantitics oi data. The need to 
obtain a mathematically tractable iV»i‘mula was a consideration in the choice of 
the fourth (quartic) power in the two-quailic profile i5(. The same mathemati- 
cal difficulty has also resulted in the existence and use of a multitude of ap- 
proximate correction formulas » 8 i. Most of these are in agreement at high 
elevation angles, but disagree at low angles, apparently because of differing 
approximations used, even when the same profile was initially assumed. For- 
tunately this discrepancy is not too serious since most satellite tracking 
measurements are taken at tlie liighei’ elevation angles. Nevertheless, it is 
possible to obtain, for any assumed profile, corrections that are both ac- 
curate and yet suitable for use on tracking data. To do so it is necessary to 
take advantage of the circumstance that the corrections ordinarily are exercised 
repeatedly using varying values of elevation angle and range but with a fixed 
atmospheric profile. Consequently, initial or ’’pre-pass” calculations which 
involve only the atmospheric conditions at the rime of the satellite pass and 
which are independent of satellite position may be lengthy without causing a 
significant percentage increase in the total computer time per satellite pass. 

For efficiency in calculation, therefore, the mathematical formulation of the 


corrections should be such that quantities functionally dependent on the 
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parameters of the atmospheric profile are separated from those dependent on 
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elevation and range. 

Such a separation can be effected by expansions in rational functions of the sine 
of the elevation angle and negative powers of the range. The coefficients in these 
expansions will depend on the atmospheric profile alone, and can be calculated 
in advance of the satellite pass. 

2. GEOMETRY AND NOTATION 

The geometry involved is shown in Figure 1. The earth is taken to be spherical 
with a radius 'a* nominally equal to 6369.95 kilometers. The tracking station is 
located at a distance r^ from the center of the earth, and at a height h^ above 
sea level. The satellite is at a distance from the center of the earth, and at 
a height h^. The radio ray path betweexi the satellite and the station is shown as 
a dotted line. The distance between the center of the earth and a given point on 
the ray path is *r'. The height of the point above sea level is 'h*, and the eleva- 
tion of the ray path at the point is S . 

The angle-of- arrival is the angle 8^ above the local horizontal at the ground 
station. The angle Ae is the elevation error, or difference between the angle- 
of-arrival and the true elevation angle E of the satellite with respect to the 
ground station. 

A E = - E (1) 
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and is the total l>ending of the radio path. The refractive index n is assumed 
to depend only on the height h alx)ve the surface of the earth. 'J1ie i adio re- 
fractivity 


N (h) io« In (h) .. 1] (2) 

will be used in normalized form, and a normalized height aix>ve the tracking 
station will be employed. 7'aking the refractive index and the refractivity at the 
tracking station to be 


c'>,) 

(3) 

(h„) 

(•i) 


respectively, and taking the effective height of the troposphere above the track- 
ing station to be 


H = (1 n;) 


j N (b) H h 

dj. 


(5) 


the normalized height above the tracking station 

X = (h - h^) H (6) 

is defined. In terms of this variable, the refractivity may be written as 


N (h) ^ Nq f (X) 


(7) 
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where the normalized profile f(x), which will be abbreviated as f, Is 


f ^ N (h^ + H xV'Nq 


which is equal to unity at x = 0 


f (0) 1 


and integrates to unity 


f d X = 1 


from the tracking station upward. In the particular case of an exponential pro-* 
file, f is equal to exp (-x). 


The radio range, or electrical distance along the ray path is designated as H , 


'• ■ J" 


(n/sin i9) d r 


the geometrical distance along the ray path is 


f'’ 

I Cl/sin 6*) dr 


and the straight-line distance or slant range is R. The range error is then the 
difference 



(13) 


r 
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AR = R - R 

t* 


The parameters 


P = v2H'rp 

(14) 

q = 10“6 Nq Tq H = 2 X 10-« N^ /'P2 

(15) 

Q - q cos2 6 ^ 

(16) 

the normalized sine of the angle of arrival 


o. ~ (1 ^p) s i n 6^ 

(17) 

the normalized sine of the elevation angle, 


0 ~ (1 'p) sin E 

(18) 

and the normalized inverse range 


= p fjj/R 

(19) 


will be used. Typical values are p = 0.05 and q = 0.25. The value of a range 
from 0 on the horizon to about 20 for a wave arriving vertically downward. 

3. INPUTS TO THE CORRECTION FORMULAS 

The desired form of the correction formulas for the elevation error A E and the 
range error AR depend on the intended use. In the first case to be considered, 
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which applies, for example, to a tracking radar, it is assumed that measured 
values of the angle- of- arrival and the radio- range are available. The 

quantities /\E and aK should then be given as functions of the variables 6*0 and 
Re- 
in the second case, the situation is somewhat different. It is assumed that the 
satellite ephemeris is known from previous tracking and orbit determination. 

In this case the true slant range R and the tioie elevation angle E will be known 
quite accurately as functions of time, and one or both of the corrections AE and 
/ R may be needed to provide accurate predictions of the angle-of- arrival and/or 
radio- range and/or range rate either for acquisition or for comparison with 
values to be measured. The latter comparisons are used to improve the satel- 
lite ephemeris by an iterative process that minimizes some weighted function of 
the differences observed. Here AE and AR need to be expressed as functions of 
E and R. 

In both cases above, the percentage difference between R and R^. is small. R^ 
is at worst about 200 meters larger than R, while R has been taken to be at least 
70 km. Consequently the distinction between these quantities may be neglected 
in their use as inputs to the error formulas. The distinction between E and 6 ^ , 
however, must be retained at lower elevation angles. 
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4. CORRECTIONS USING KNOWN ANGLE OF ARRIVAL 


4.1 Formula for Elevation Error 

A procedure for determining the elevation error has been given by Bean and 
Thayer [9], The bending t is first calculated by integration, and the angle .• is 
next calculated from the geometry of Figure 1 using the known value of r . The 
elevation error is then given by 


A E = T - 6 (20) 

Using the notation of Section 2, the usual expression for the bending is [2] 

T :r 10“® Nq cos I ( a) 'p 
where the bendii^ integral is defined as 


I(a)= r Hv (22) 

Jij \/x + - Q (1 - f) 

in which it has been assumed that the satellite height hi is great enough to per- 
mit the upper limit to be extended to infinity. 

The equation giving o in terms of r is derived in Appt ^ix 1. Using (21) and 
(1-4) in (20), the formula for the elevation error is 

A E = 10"3 Nq Cl ^P) cos [I (a) - p L (a)] mrad (24) 

with 


L ( a) - 1 - a I ( a) 


+ iq (a.) 

4 


(25) 
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4.2 A Digression on the Approximations Used 

The quantities 10“'’ ' 300 x 10“^ and p^ 0.0025 are neglected compared to 

1. However the quantity q ' 0.25 cannot be so neglected, even though it contains 
the surface refractivlty as a factor. 

The quantity Q that appears in the radical in Equation (22) is not independent of 
satellite position because of its dependence on . This would lead to complica- 
tion of the formulas to be derived. Fortunately Q may be replaced by q with 
negligible error, since the term neglected thereby is small compared to the 
square of - 

q sin2 (1 f) << a2 - (1 p2) sin^ 0 ^^ 

It will also be found that the same approximation canS and in some cases 
should^, be made elsewhere in the formulas to be derived. 

4.3 Expansions of the Bending Integral 

In order to use the formula (24) for the elevation error a rapid method for cal- 
culating the bending integral. Equation (22), is needed. Since a is as large as 
20 at liigh elevation angles, it is natural to expand (22) asymptotically^ in powers 

1. In .’•hi? asymptotic expansions that follows, when Q=q - qp^a^is set equal to q, the term 
neglected is small compared to the preceeding term in the series. 

2. If the approximation is used in calculating I (o^ then it should be used in calculating 
polynomials involving I(o0 (for example, in eq 25) in order to obtain the correct asymptotic 
expansion of the polynomial. 

3. The expansion is asymptotic it q = 0 and f = exp (— x). The nature of the expansion was 
not investigated for other cases, Pnd the development proceeds on a formal basis. 
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of l/a . This can be accomplished either by a formal binomial expansion of the 
radical in Equation (22) after the square of a has been factored out, or, alter- 
natively, by repeated integration by parts of the numerator of the integrand. In 
either case there results (see Appendix 2) 


I (a) - (l a) - (l/a)3 f (l/a)5 


(26) 


where 


I. -i 


(27) 


and 


I, = (3/4) 


[( 


X f d X - q ( 1 


j fa, 

0 


X + (1/6) 


(28) 


Use of the first term of (26) in (21) results in the familiar 


T - 10”^ Njj cot tnrad (29) 

which holds at high elevation angles. 


Equation (26) suffers from the usual property of an asymptotic expansion - it is 
not useful at small values of the argument a. There is, it happens, a procedure 
[10] for converting a divergent series such as Equation (26) into a continued frac- 
tion expansion that converges, in this case, for all a > 0. The expansion diverges 
at a = 0, however, and converges only slowly when a is near zero. Rather than 
apply the procedure directly, therefore, the integral I( a) is expanded, with 


small, as 
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(30) 


where (Appendbc 3) 


p 

i„ 1(0) 2 V X q a n - -- .lx (31) 

(1 i. q I')* 

and 

i J I' (0) -2 r (U) il li q f' (0)1 (32) 

Noting that Equation (26) approximates I( -) when is large, and that Equation (30) 
approximates I( < ) when / is small, the approximation of I( .) over the entire 
range of a is accomplished uy means of a ratio of polynomials in ■ . The coef- 
ficients of the polynomials are chosen in such a way that the expansion of the 
ratio in inverse powers of .. agrees with the leading terms of Equation (26) on 
the one hand, and its expansion in ascending powers of t agrees with the leading 
terms of Equation (30) on the other hand. Tins method of approximation insures 
accuracy if a is either large or small. Accuracy with intermediate values is 
obtained by the inclusion of a sufficient number of terms from each series ex- 
pansion, The number of terms used here — three from Equation (26) and two 
from Equation (30) — is not necessarily optimum, but worked out well when the 
method was applied to an exponential profile. It is evident that the method re- 
quires a certain degree of smoothness in the profile if an accurate approxima- 
tion is to be obtained with this number of terms. 
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4.4 Form of the Approximation 

Consider a rational function of a, f( x; F, , F^ , f^, f^ ). which depends on four 
parameters Fj , Fj , f ^ and f ^ , and which is expressed in the form of a continued 
fraction 


FC: F,. F,; f„, f,) = 


a + 


where the intermediate constants J[i,f 2 ,^ 3 , and_f^ are calculated from the set of 
parameters Fj , F^. f o» and fj using in sequence 


f. = F. 


i. --(F,/f,)-f, 




i*= f, f, 

— 4 0 —1 —3 —2 


On clearing the denominator of Equation (33) of xractionsy and expanding the re- 
sulting fractional form by long division in descending powers of a 
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F (a; Fj, Fji fj) r; (1 /a) - Fj (1 0^ * F, (1 .) 5 -• constnnt •(! /v* • •(38) 


If the long division is carried out using ascending powers oi 


F (a; Fj , F^ : fy fj) “ - fj > ^ const ant 


Thus the function F(a; Fj , Fj , fp, f j) is wGll-suited to approximate I( ), Equa- 
tion (21), provided that the parameters Fj, F 2 , f,,, and fj are chosen as I, , l^, 
io, and ij respectively, i.e., 1(a) is approximately equal to F( >; Ij, ij ). 

Note that these latter parameters. Equations (27), (28), (31), and (32), depend 
only on the refractivity through the parameter q and the profile f(x). Their 
calculation is independent of satellite position, and numerical integration prior 
to each satellite pass need not be ruled out. If f is a given model profile, more- 
over, the integrals can be evaluated analytically if the functional form of f(x) 
permits, or, otherwise, may be evaluated numerically and curve or surface 
fitted empirically. 


4.5 Formula for Range Error 

The range error Equation (13) may be written as the sum of the difference be- 
t^veen the electrical and geometric distances along the ray path and of the differ- 
ence between the geometric distance along the ray path and the slant range 


A R = (R^ - p^) . (R^ - R) 



The first term ii* Equation (40), the difference along the ray path, is, after 
division by r^ 


! 1 


H 


(Re - '^0 U 'r,,) j* ( 10 "^^ N (h) /sin (j d r ^ i 10 "<^ N^,p J (a) ( 41 ) 

Jo 


where 


J (a) 


[■ 

Jo V 


X 4 - q ( 1 -- f ) 


d X 


(42) 


9\ 
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While the geometrical difference (R^ - R) between the ray and straight-line paths 
is smaller than the electrical difference (R^ " g) along the ray path, it is not 
entirely negligible. In Appendix 4 an expression for this difference is derived. 
Substituting Equations (4-8) and (41) into (40), the expression tor the range 
error is 


AR/To =^10-«N^P [M(a) -lpOLJ(a)] 


(43) 


where 


M ( a) = J ( a ) + q 


I (a) - 1 K (a) - i a 12 (a) + i q 12 (a) 


(44) 


with 


- .7,^ ^ f' 

0 Vx + - q (1 - 


d X 


(45) 


f) 



The expansions of J( a) and K( a) for large and for small values of a have been 
given in Appendices 2 and 3. Thus it is possible to calculate these integrals 
in the same way as was done for I( a) by using Equation (33) with suitable choices 
for the parameters Fj , Fj, f^ , and f^ and then to compute M( a) by substituting 
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the calculated values of I( a)» J( a) and K( i) into Equation (44). It is more efficient 
to calculate M(^) directly however. Thus, substituting Equations (26), (2-3), and 
(2-6) into (44) 


M 


( ) - (l/o) - Mj (l/a)3 + Mj (l/a)5 - 


where 


M, = L 
* 2 


O) 

0 


X f d X - q f 1 


-if 

•'o 


d X 




jij x*fdx-q^lf j|^ xfdx-i J xf^dx 


(46) 


(47) 


(48) 


Substituting Equations (30), (3-7) and (3-11) into (44) 


M ( a) = hIq - iHj a -f . • . 


(49) 


m. 


ni. 


)o +Vo + (1^12) q2 ij - iqk^ 


= jl 


(50) 


(51) 


M(a) is now calculated using Equations (33-37); i.e., is given by F(a; M^, M^, 

m^.mj). 


A formula for the range- rate correction can be derived (Appendix 5) and, 
presumably, modelled as above, but this has not been done. In practice (with 
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the Goddard Range and Range Rate System) the range- rate correction can be 
obtained by dividing the difference between successive range corrections by the 
corresponding time interval. 

5. CORRECTION'S USING KNOWN ELEVATION ANGLE 

The formulas for the corrections AE and dR as functions of d and R rather than 
as functions of « are more difficult to obtain Reichly 14] in an excellent treat- 
ment of the problem, has used a perturbation method to obtain the first two 
terms of the expansion of these corrections in powers of the surface rofrac- 
tivity. This, in the notation and with the approximations used here, is equiva- 
lent to an expansion in powers of the parameter q which may be as large as 
0.64 at Nq = 450. Such an expansion is accurate at large values of the elevation 
angles as may be determined by an examination of the asymptotic expansions, 
in which the higher powers of q appear only in the higher order terms. More 
terms are needed, however, if formulas numerically accurate at small values of 
the elevation angle are to be obtained by this method. Such formulas can be ob- 
tained by expansion, not in powers of q, but in powers of p, 6 and l/d, and by 
the use of the method of the preceding paragraphs. 


5.1 Elevation Correction 

It is assumed that R is large enough to permit A E to be calculated from the 
first two terms of its expaiision in powers of the normalized inverse range p . 
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( 58 ) 



the parameters and are determined by substituting Equation (58) into 
Equation (56) with I in (56) given by Equation (26), expanding in powers of 1/ , 
and equating coefficients of like powers. The results are 


k ! Jfh 


xfdx+q(i.4l f 

2 J, 


d X 4. (1/6) 


At small values of 


U (/O = Uq - Uj /5 + u. 


/32 _ . . . 


To obtain u« , 6 is set equal to zero in Equation (56) to obtain 


% = % 


which can be solved iteratively. The parameters u. and u- are obtained by 


differentiating Equation (56) 


whence 








1 - Iq I' (^Iq U 


Similarly, a second differentiation results in 






Using U|i U 2 , Ug, and in Equation (52) is approximated by 



F r/3; . Uj ; Uq, u^). 


To calculate V(/3) Equation (51), Equation (63) is solved for I' (6 + 1/2 q u^) and 
substituted into Equation (57): 

V 1 - /3 U ( 0 ) - 1 q fl2 (/3) 1 + i q U’ ( 0 ) ( 66 ) 

Here U'(/3) may be calculated using , Uj, u^, and u^ as shown in Appendix 6. 


5.2 Range Correction 

Substitution of Equation (54) into (43) and expansion in powers of ,0 gives (using 


Equation 5-1) 


^ '^ ''■0 = 5 10-6 ^ M ^,5 + i q t i(J p L2 i q uj 


which becomes 


* ^ N„ p |w (ij) t 1 Q p 1 - ,S U (S) - 1 q (/?) 'I 

where, using the asymptotic expansion of U in M ( /5 + 1/2 q U) as given by 


Equation (46) 


( 0 ) - ( 1 / 6 ) - Wj (1/^3 + (l//3)5 - . 




X f d X + 


K 



x2 f d X f (q/6) 1 I 3 


f 2 d X I f f f ^ d X (*^1) 


and where, setting 0 ~ C 


W (/•) “ Wq - Wj 


P f • • • 




- u; 


Wj = 2 1 




6. EXAMPLE USING AN EXPONENTIAL PROFILE 

If an exponential refractivity profile is assumed, the normalized profile f be~ 


comes 


f (x) e' 


The effective height H is equal to the reciprocal of the decay constant of the 
profile and will be estimated from the refractivity at the station using the em- 
pirical formula [10] 


I'H = In 


M -T o/x 0-005577No 

Nq - 7.32 e ® 

6.1 Known Arrival- Angle 

Most of the required integration, that in Equations (28), (47) and (48), may be 
performed directly. Equation (3-5) or (31) and Equation (3-12) however, must 
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be handled numerically. The simplest procedure is to evaluate these integrals 



numerically [12] at, perhaps, 20 different values of q spread over the expected 
range 0.1 < q <0.7. The values obtained may then be approximated by a poly- 
nomial using, since it is simple to apply, a least squares fit. One finds, for 
example, that a cubic polynomial in q approximates i^ to better than 1/4 percent 
over the range 0 to 0.7. As it happened, however, an exponential expression 
gave a better approximation with fewer empirical constants: 

Iq v'^d - 0.9206 qr0 '*^6® ± 0.04% 0<q<0.7 (77) 

The integral was also closely approximated by a similar exponential 
expression. 

The equations for the corrections are collected in Appendix 7. 

6.2 Known Elevation Angle 

Here, also, the coefficients of the asymptotic expansions, Equations (60), (70), 
and (71), may be integrated. Equation (62) was solved numerically for at a 
number of values of q by Wegstein*s method [13], using repeated numerical 
integration. The calculated values of u^ were fitted by 

Ujj = \/;r(l + 1.4844 q) - 0.39 144 + q. 02%0 < q < 0,7 (78) 

The coefficients Uj and u ^ were also obtained numerically using the same cal- 
culated values of Uq in Equations (64) and (65) with I' and I" given by Equations 
(3-3) and (3-4). 
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The equations are collected in Appendix 8. 


6.3 Numerical Examples and Comparisons 


Let the refractivity at the tracking station be Nq =313. If the tracking station 
is at sea level, r^^ = 6369.95 kilometers. Performing the pre-pass calculations 
of Appendix 7, H = 6.951 km, p = 0.04672, q = 0.2868, and the equation for the 
elevation error is 


with 


E = 0.313 cos 6 ^ [i - (6369.95/R) L] mrad 


0.00093^8 


sin + 


0.002117 


sin +. 


0.006054 


sin +. 


sin + 0. 1163 


L = 1 - i sin + 0.0001565 i^ 


The equation for the range correction becomes 


A R = 0.002176 [m - (913.5/R) cos^ $^] 
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1 

m 

0.0008565 

+ 

0.002173 

sin + (®3) 

0.006082 

sin + 

sin 6'^ +0.1157 

Using twelve values of from 0 to 900 mllliradians and two different values of 
the range R (in kilometers) at each value of (these ranges are for satellites 
at heights with respect to the tracking station of 70 km and 475 km), the cor- 
rections in Table 1 were calculated. The corrections calculated using a double- 
precision ray trace program [14] are also listed for comparison,® 

The largest difference is 0.3 percent. If desired, a final empirical adjustment 
of the coefficients can be made to reduce this error (Appendix 9), In Table 2 
the corrections calculated from Appendix 8, using the values of elevation angle 
and range listed, are compared with those obtained from the ray-trace program. 
Since in this case the corrections do not hold for negative elevation angles, all 
but the smallest of these have been omitted from the table. The maximum error 
in Table 2, 0.9 percent, is larger than that found in Table 1, apparently because 
only two terms were retained in Equation (52). 

6. These values will be found to differ somewhat from those in the CRPL roy-troce tables [15] 
which used single precision and slightly different values of rQ. 
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Table 1 

Calculated Corrections, Arrival Angle Known 
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Table 2 

Calculated Corrections, Elevation An^e 
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Similar checks were made at 200 and 450 for both kinds of correc- 
tions with equivalent results. 

The computer time required to perform the calculations is also shown in the 
tables. The pre-pass calculations in Table 1 took two sixtieths of a second 
(IBM 360/95). The two sixtieths of a second used to compute the twenty 
four pairs of corrections in Table 1 rtjpresents a rate of about 700 per sacond. 
Programs designed for operational use should show improvements over these 
times. 


7. SUMMARY 

A method is given for deriving atmospheric elevation-error and range-error 
correction equations in a form suitable for use in the processing of satellite 
tracking data. The method can be applied to any sufficiently-smooth spherically- 
symmetric model of the atmospheric refractivity. The method was tested by 
application to an exponential profile and comparison of the resulting corrections 
with those given by a double-precision ray trace program. The results were 
in agreement to better than one percent over the entire range of elevation 
angle (0 - 90°) and better than 0.3 percent over most of the range. 

An inspection of the coefficients in the correction equations shows that the cor- 

i! 

rections at high elevations do not depend on the final detail of the refractivity 
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proflI< » li(h), but rather on the value of the surface refractivlty and, in apparent 
orde:* :*li njwrtance, the height integrals j Ndh, J hNdh, and J N^dh, etc. The 

deriv? N'(0) of the refractivity at the tracking station is also significant be- 
cause it v'l^ighs heavily in the determinatit a of both coefficients of the small- 
angle e xpansion of the bending integral. Where improved accuracy is required, 
emph(.;4i:? sliruld be placed on studies leading to better estimates for these sig^ 
nificam r> entities. 
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APPENDIX 1 


EXPRESSION FOR IN TERMS OF / 


/ 



\ 


! 


The equation for the angle follows from the geometry of Figure 2, and from 
Snell's law for a spherically stratified medium. In Figure 2 r^ , r^^, and R are 
all projected or a line peri>endicular to the tangent to the ray path at the satel- 
lite. It follows that 


cos fy cos - t) f R sin c (1"1) 

Snell's law for a spherically stratified medium is 

n r cos H - Hq cos (1“2) 

At the high altitude of the satellite n is equal to unity, r = r^ , and e - Com- 
bining Equations (1-1) and (1-2) gives 

sin b “ (Tq-'R) [n^ cos 6 ^ - cos ( 8 ^ - t)] (1"3) 

Approximating sin b by b, sin t by r, and cos r by 1 - 1/2 r^, 

S = (To /R) 


i’ 


1 

0 2 


cos 


- r s i n 


0 


(1-4) 
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Figure 2- Projection of R, Tq, and r. 






APPENDIX 2 


ASYMPTOTIC EXPANSION OF I, J, AND K 


Setting Q - q in Equation (22), factoring a out of the radical, and expanding 
the latter by the binomial theorem 


I (a) - a /a) 


*>-0 C 

-f'ji-i 
Jo 1 2 


[x - q (1 - f)l (l/a)2 


(2-1) 


+ (3/8) [x - q (1 - f)]2 (1/a)^ - . . d 


The terms in (2-1) that do not contain x as a factor may be integrated directly. 
Those that do contain x can be integrated by parts on the factor containing f* and 
powers of f. The result of these integrations is given in Equations (26), (27), and 
(28). (The same result can be obtained by repeatedly integrating the numerators 
starting with Equation (22), using the definition 


f(-i) 


f d X 


and assuming that f (“D vanishes with sufficient rapidity at large values of x). 


The integral J(a) Equation (42) may be expanded in the same way with the result 


J (a) - (\/a) - (l/a)3 + (ij/a)^ - . . . 


(2-3) 


34 




J - - 


r 


X f d X - q 


(■■ I ""jj 


J, = (3/4) 


U 


x^fdx-qfj xfdx- | xf^dx] 

'o ‘0 


4 q- 


'i. f 

(2 h 


f2 d X 4! 

2 


d X 


(2-4) 


(2-5) 


Similarly Equation (45) yields 

K (a) ~ (1 / a) - Kj (1/a)^ 4 (1/a)® 


(2-6) 


‘.4f 

*/o 


f2 d X - (1/6) q 


K, = (3/4) 


(2-7) 


r r°° 

-a. r'" 

1 

' j X f2 d X - q 

f2 d X - (2/3) f3 d X 

+ (1/12) q’l 

I^Jo 

Jo Jo 
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APPENDIX 3 


INTEGRALS FOR I, J AND K, AND 
EXPANSIONS FOR SIAALL VALUES OF a 

The integral Equation (22) can be integrated by parts. Multiplying the numer- 
ator and denominator of the integrand by 1 + qf’ and integrating the radical by 
parts 


I (a) = a . 2 PvT 

1 + q f' (0) I 


+ a- 


2 - q a - f) 


■ 


(1 +q 


d X ( 3 - 1 ) 


The first term may be absorbed into the integral 


1 (a) ~ 2 I [vx + - q (1 _ f) _ a] 

Jo (l+qf')2 


f" 


d X 


(3-2) 


Differentiating Equation (3-2) 


I' (a) = - 2 


r 

Jo 


1 - 


a 


v/x + - q (1 - f) 


f- 


(1 +q f ')2 


d X ( 3 “ 3 ) 


Differentiating Equation (3-3) 


I" (a) . 2 


X - c 
• 0 [x + — 


q (1 - O 


f" 


d X (3“4) 


q (1 - f )]3 '2 (1 ^ q f ')2 


I 
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Equation (31) for 1(0) is obtained by sitting a equal to zero in Equation (3-2). 
The integral 



I (0) = f r V d X (3-5) 

Jo \ X - q (1 - f) 

obtained directly from Equation (22) can also be used, but more care must then 
be exercised since the integrand diverges at x = 0. 


Equation (32) follows from Equation (3-3) by settii^ a = 0 and integrating. In 
the same way, from Equation (42) 



37 




and netting 


K (a) = k,j - kj a f . . . 


k 


0 


- a 



“ q (1 - h 


f f** ^ f q f 

(1 > q f')2 


k, = -4 f' (0)/[l + q f' (0)1 


2 i, 



APPENDIX 4 



APPROXIMATION FOR THE GEOMETRICAL DIFFERENCE 


Using Ecyuation (1-2), R^ equation (12) may be written as 

R - p ♦ »o ^0 ^0 “ »o »o ^ 0 J ^ (4^1) 

Jfo v'r2 - (nj rj cos2 t\)/n^ 

The left-hand term in the integrand can be integrated. On setting n( r^) = 1 and 
making use of the geometrical relationship 


sin R cos (5 + Tq sin - t) 


(4-2) 


which follows from Figure 2, and using the exact integral for r [9] 

to. 


r - - COS I . d r (4“3) 

Jq n2 V - n2 r* cos^ 


there results 


(R^ - R)/rQ = - (R/Tq) (1 - cos S) -(1 - cos r) sin 6^ 


(4-4) 


+ (r - s in r) cos 6 'q + cos^ 


p I 

J 1 

•'n a/ 


"o 


d r 


1/ n^ r? 0052 P ) 

n3 y r2 . ° ° ° 


n“ 
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On setting 


cos 0 
s i n r 


cos r 



(4-5) 


(4-6) 


(4-7) 


and making the usual approximations in the integrals, there results 



R) 


i 10-« p O 



iK-i 

2 2 


12 f q 13/12 


- (Tq p/2 R) 


1 



(4-8) 


APPENDIX 5 



RANGE-RATE CORRECTION 

The range- rate correction is the difference between n(r j) times the projection 
of the satellite velocity on the ray path at the satellite and its projection on the 
straight line joining the satellite and the ground station. 

It may be found by differentiating Equation (44) with respect to time. Using 
primes to denote differentiation with respect to a , and the dot notation to denote 
differentiation with respect to time, and making use of the relation (obtained by 
performing the indicated differentiation and then integrating by parts) 


i' 





appends; g 


CALCULATION OF U'( / ) 


From Equations (59) and (GO) 


U' - ^ 3 U, (1 - * 5 (\n^ • 

U’ - [Uj - 2 Uj /3 t . . .J 

The calculation of U* is accompliBhed by means of the form 

G (0) \ 

o2 ^ 

M2 

1 f 


(0-J) 






where 


= Gj (G“ 4 ) 

S2 = “ Ki 

£ 3 = S/(gigo- 

S4= £§Rj g,/g2 

whence 

G r il/3)^ - Gj (\/0)^ 4. (l//3)6 . const. //37 . . . ^;-h) 
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and 


G - gy - /: f constatit i . 

Setting Cr , = 3Uj , G2 = SU^, g„ = u j, and g ^ = 2U2 
equal to minus G. 


(6-9) 

IT' will be approximately 


APPENDIX 7 


CORRECTION EQUATIONS FOR AN EXPO^'IENTIAL 
PROFILE, ARRIVAL ANGLE KNOWN 


H r l/ln 


\ - 7.32 eO-OOS577 Nq 




q = 10-6 


= v^(l - 0.9206 q)-o-4-»68 


ij r 2/(1 -q) 

^ 2 \ 2 7 

I 2 = 0.75 [1 - 0.75 q +(l/6)q2] 


\ = v/2T (1 - 0.9408 .irO'-^^sg 






1 _ . 


"'l " 2 1 +;jq /(I -q) 




.. . 0.75 1 - s, . !1, 


i = F (sin Ij. i()/p. ij./p^' 

L = i - i sin + -10-6 N i2 
6 2 ° 
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m r i ij 


/E 0.001 cos ^ (i - L R) 

ni F(sin ^ ni^^/p , ni^/p^) 

- i 10"fi cos2 - , L 2 H) km 

2 0 0 j 

The function i above is equal to I( a)/p, and the intermediate constant, Equations 
(34-37), associated with i can be obtained from those of I( a) by multiplying the 
first three by p^ and the last by p. 

The programs used to perform the calculations were written in A PL language 
[16], and are shown in Figure 3. 


.R - 10‘® N„ H 



VPR EPASS CA LCVLATIOPSimv 

V PREPASSCALCULATIOESl 

[ 1] HH^r%NNO . 32x*0 . 005 577xA’;V0 

[2] PSQVARED<-2^HIHR0 

[3] P*-PSQUARSD*Q .5 

[4] 

[5] F0-^70+-( (oi )*0 . 5)x( 1-0. 92 06 xf?) *”0.44 68 
C6j F1^I1^2\1-Q 

[7] Ff . 5x1-0 . 5x(j 

[8] FF2^II2-^0 . 75 xl t^x”o . 75+C f 6 

[9] CALCULATEl 

[10] 11^PEQUARED>^IX 

[11] 12^FSQ7ARED^12 

[12] l^i<-PRQUARED^gZ 

[13] X4^PxF4 

[14] A'0*h( (02 )*0. 5 )x( 1-0. 94 08xQ)*“0 . 47 5 9 

[15] F0^A/0^( J0xl + ^xi+(-^xl0xJ0il2 )-QxA04 2 

[16] F1*^A'1-<-2x( l + (^xJ0xJ0*4)il-g 

[17] FFl^i^/Ail^O. 5x1 -0.75 x(? 

[18] FF2^MM2*-0 . 7 5 x 1 + (,ix ( “2 5 t 24 ) + 1 1 x^ f 3 6 

[19] CALCULATEF 

[20] M1^PSQUARED>^F1 

[21] M2^PSQVARED>^F2 

[22] m^PSQUARED^^gZ 

[23] M^4^PxgH 

V 


V CORRECT TORS 1 CD] V 

V CORRECTIONS! 

[1] JIF^lo7FF7/l Ox 0.001 

[2] C0S^2oTHETAQxQ ,00! 

[3] J-«-vFJA + I1^JJA’+I2tFJF+I3 ^5I//+I4 

[4] LL*-! -Ix;^Iu-o . b xi/-;“6x,v;'0xj 

[5 J DLE<-(J . OQlxNNQxCOSxJ -RQxLL^RR 

[ 6 ] M^-iSIN+H!TSIN+tl2 iSIN+I^Z ♦5IA+i/4 

[ 7 ] DRR^!E''ZxNN0xHHxM-0 , Sx!E' exNlJQxROxRQxCOSxCOSxLLxLL^RRxHH 

V 

VCALCULATEFUJ2V 

V CALCULATEF 

[1] Fl^FFl 

[2] F2^{FL 2rFl) -F! 

[3 J F3-<-F2 v(; uxF0xFlxi+£lfF2)-l+FlxFl 

[4] F4-i-F0xF1xF3^F2 

V 


Figure 3. APL Programs for Corrections, Arrival Angle Known. 
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APPENDS. 8 


CORRECTION EQUATIONS FOR AN EXPONENTIAL 
PROFILE, ELEVATION ANGLE KNOWN 

H, p, and q are calculated as in Appendbc 7. 

2 (1 4 1.482 q )-o 38 26 

2 v“(l f 1.71 q)0 1 


Uq ^ \iv (1 + L4844 qr^ 



(1 + 1.6454 q)-o 583 




r 


u r F (sin E; Uj. U^; ^ij/p^) 


u' = - G (sin E; 3 p^ U^. 5 p^ U^; n^ ''p^. 2 u^'p^) 

V ^1 - u sin E - 1 10“® Ng (1 » 10“® u') 

iE - 0.001 Nq cos E (u - V Tq /R) rnrad 


F (sin e; p^ W^, p** W^: /p. ^p^) 


A R - 10“6 Nq H w - i 10^6 Nq ~ cos2 E fl - u sin E - i 10'® Nq j km 


The corresponding A PL program is shown in Figure 4. 


V .li-iV/'ViJ J I U J •* 


[1 J 

HH'*- v»PPy+pp'U-7 . 3 

3 »■ *■ sj . 0 0 5 5 7 7 Xj 0 

[2 J 

ic L' w' ^ •’ » * '* * '■■ *■ ' '■ •' *■ ’* 


[3] 

F^-PCmjAUi:i)*Q . b 


[4] 

OxpPu+PUvPE 


[5] 

IT'?/-*- 3x( 1 + 1.482x 

p ) + 0.382b 

Lb J 

DEE^l -0 . b^C'i^S FU 


[7] 

JPPJ*-2*( r 01 ) *0 . 5 

) X ( 1 + 1 . 7 1 X 'i ) ★ 0 . 1 

[8] 

F(}4~ij\j*-i (01 )*u . 5 ) 

X ( 1 + 1 . 4 8 4 4 X , ~ 0 . 3 9 1 4 4 

C9J 

7'i>^Pi-^--jpp ^Pi.■p 


Lie] 

FFi*‘UUl-<-0 . 5 >'• 1 + 0 . 


LllJ 

PPO'** UU2 •*~U .75*1 +£,, 

X ( 7 V i 2 ) + C V & 

[12] 

C/lIiPPE/1 P/;. P 


[13] 

P 1 PJ Q UA. P/.'tPx^l 


[14] 

PEP*P2 


[15] 

P3*-PPaP/lPEP^P3 


[16] 

k4*-p*P4 


[17] 

u3*pp--u*-pl 


[18] 

o 1*-P1'1 /-PP T P 

I'M* A 

[19] 

Or (/ 1 •*“PPPl **3 * tP7 1 


[ 20 ] 

2 ■*•[! U F2 *- J X tj U 2 


[ 2 1 ] 

C/UPPP4V'Ef2 


[22] 

ypi*-po., J rp P'oi 


[23] 

P P 3 F' J w' .1 j‘ : ■ 2 * X 0 2 


[24] 

2FA-*- Fi-C U Ai\ ED *P3 


[ 25 ] 

EP'4x-P>-o4 



[ ,<^uj KU*-( (.O^)*0.5)x(l + l.tj454x<;)*‘‘o. ^33 
[ 27 J F 0 <♦'}:* 0-»-( U0*-ltQ-^l-Q>^U0^U0ii6)-0.b^Q*-KU 
L28J F1+7/1-2X1 -0 . 25xi?xc/0xL/0 
C2yJ FFl^k'lv’l-^U. 5*1+0. 28 x£? 

L 3 0 J f F 2"'' i. k 2 0 . 7 5 * 1 + v.j X ( 7 t 4 ) + t,j t 1 8 

C31] CALCULATFF 
[3 2] k:l-^FFi2y/4A'FF*£i 

[3 3J W2*‘PFQuAFll''<F2. 

C34j WA-^PFQUAFFD^-F'i 
[3 5] 

V 


LU 
[ 2 ] 
[3 J 
[4] 
[5 J 
[b] 

[7] 

[ 8 ] 


[7 ] 
[ 2 ] 
C3] 
i 4] 


VCOFFifCT ' UFF 
C0r:RECTI0nF2 
FTFF‘^loFF*0 . uOl 
C’0FF-^2OFFx 0 .001 

U*-iSIFFi-ill iFlFE + U2 iFIFE + lJ3 iSIHE+Uk 
l/F— J ( F/F2'xFJ*VF ) +ypi + 1 +^P2 i 4fP3 +J/F"*F’ 

P^PFP*-1 -pxFJF;:: + 0. SxlF"6x’’P0*P 

PEE<-0 . 001 xEE0*C(7EExJ-P/lPEPx( 1 + lE”bxEPO *EP) xPO+PP 
W*-iSinE + Wl rSIEE^W2 iSXEF*l£3 +5JPP + E4 
EE’P«-lE~6*A^V0x////xk'+0 = 3 * 1 E” 5 * /.‘PO x ( iPH^EF 


P 4 + . 


) x(P0xE0EExP4 


VC/lEEEE/lTEk’lG]V 
V CALCULATE 2 
kl '♦'ucJ 1 

k2*-(EJ2?kU-2l 

C3*‘k2 1 l + olx(70 

k4'*"k3x^3x^lxiyi iG2 


Figure 4. APL Programs for Correcticns, Elevation Angle Known 


APPENDIX 9 


IMPROVED FIT 

The correction equations given in Appendix 7 can easily be modified to give a 
closer fit to ray- trace calculations. The method chosen was to modify the co- 
efficients f,, f , and f using a least squares adjustment that brought 

— 2 —'3 —4 

F( 1/2, 3/4; v /ft 2), /which is the approximation for 1( a) where q ~ 0) into 

2 

better agreement with its theoretical value e"' erfc a . The factors by 
which f^, and f^ are multiplied are 1.08885, 1.320903, and 1.21313 respec- 
tively. The resulting corrections are shown in Table 3. 




Table 3 

Calculated Corrections, F Modified 
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